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ABSTRACT 

We present a new determination of the dust content and near-ultraviolet/optical ex- 
tinction curves associated with a sample of ^8300 strong, W A2796 > 1A, Mgll ab- 
sorbers, redshifts 0.4<z<2.2, identified in Sloan Digital Sky survey (SDSS) spectra 
of quasars in the DR6 release. Taking into account the selection effects that re- 
sult from dust extinction, including the reduction in the signal-to-noise ratio of an 
absorber appearing in a reddened quasar spectrum, we find a stronger dependence 
of E(B — V) on absorber W A2T96 than in other published work. The dependence 
of the median reddening on W^ 2796 can be reproduced by the power-law model 

E(B - V) = ((8.0 ±3.0) x 10- 4 ) x (W ) (3 - 48±0 - 3) for 1.0<W <5.0. Observed Mgll 
samples, derived from flux-limited quasar surveys, are shown to suffer from significant 
incompleteness at the level of 24 ± 4 per cent for absorbers with Wo > 1 A and 34 
± 2 per cent for absorbers with Wo > 2 A. Direct determination of the shape of the 
near- ultraviolet extinction curves, using high signal-to- noise ratio composites, for ab- 
sorbers as a function of E(B — V) show evidence for systematic changes in the form 
of the extinction curves. At low E(B — V) (<0.05), the extinction curve is feature- 
less and well represented by a Small Magellanic Cloud (SMC)-like extinction curve. 
For intermediate E(B — V)s (<0.2), approximately a third of Mgll absorbers show 
evidence for a 2175 A feature and an extinction curve similar to that of the Large 
Magellanic Cloud (LMC). For the small number of high E(B — V) (>0.3) absorbers, 
the majority of which exhibit strong CaIlAA3935,3970 absorption, there is evidence 
for the presence of a 2175 A feature as strong as that found in the Milky Way (MW). 
Near-infrared photometry for six of the systems indicates that the rest-frame optical 
portion of the extinction curve for these high- E(B — V), and likely very high column 
density, systems is significantly greyer than the SMC, LMC or MW extinction curves. 
Application of the new results on the dust content of strong Mg II absorbers shows that 
dusty absorbers can account for a significant proportion, up to a factor of two, of the 
observed overdensity of absorbers seen towards Gamma- Ray Burst (GRB) sightlines, 
compared to sightlines towards quasars in flux-limited samples. 

Key words: dust, extinction - galaxies: ISM - quasars: absorption lines - gamma-ray 
burst: general. 



1 INTRODUCTION 

Since the pioneering work of Bergeron (1986), Lanzetta & 
Bowen (1990, 1992) and Steidel et al. (1994) the presence 
of Mgll AA2796,2803 absorbers in 'haloes' extending to dis- 
tances of ~50-100 /i -1 kpc about luminous galaxies has been 
well established. The more recent recognition of the impor- 
tance of outflow, infall and feedback processes in general 
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for our knowledge of galaxy evolution has reinvigorated at- 
tempts to understand the processes responsible for the exis- 
tence of extended gaseous haloes associated with luminous 
galaxies. 

Strong Mgll absorbers, Mgll A2796 rest- equivalent 
width (EW) >0.5 A, seen in the spectra of background 
quasars reveal the presence of relatively cool, T^IO 4 , ionised 
gas with neutral hydrogen column densities of ~10 18 - 
10 22 cm -2 . The availability of very large samples of interme- 
diate resolution, moderate signal-to- noise ratio (S/N) quasar 
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spectra from the Sloan Digital Sky Survey (SDSS) (York 
et al. 2000), has allowed the compilation of extended sam- 
ples of Mgll absorbers (e.g. Nestor et al. 2005; Prochter 
et al. 2006; Quider et al. 2011) with consequent improve- 
ments in our knowledge of the statistical properties of the 
absorbers as a function of key physical parameters. 

The homogeneity of the SDSS spectroscopy coupled 
with the quantitative nature of the quasar target selection 
(Richards et al. 2002) has also enabled progress to be made 
in quantifying the dust content of absorber populations at 
cosmological distances. Wild & Hewett (2005a) first demon- 
strated that the Call AA3935,3970 absorber population, a 
subset of very strong Mgll absorbers, showed evidence for 
the presence of dust with E(B— V) values as high as 0.2 mag- 
nitudes. The strong Mgll absorber population as a whole 
shows much lower dust content, with mean E(B — V)~0.02 
(Wild et al. 2006; York et al. 2006). More recently, Menard 
et al. (2008) quantified the dependence of the dust content 
of Mg II absorbers as a function of both absorber EW and 
redshift. 

Interest in the properties of dust associated with Mg II 
absorbers has grown with the explicit linking of the absorber 
population to the star formation history of the galaxy popu- 
lation (Menard et al. 2009) and the association between very 
strong absorbers and outflows (Bouche et al. 2007; Nestor 
et al. 2010). Dust has also been suggested as a possible ex- 
planation for the puzzling discrepancy between the increased 
incidence of Mg II absorbers detected in Gamma- Ray Burst 
(GRB) compared to quasar sightlines (Prochter et al. 2006; 
Porciani et al. 2007). 

In the local Universe the principal observable difference 
between dust extinction curves is the presence of the 2175 A 
feature, the strength of which is inversely correlated with 
the increasing rise of dust extinction into the ultra-violet. 
The 2175 A feature is strong within the Milky Way (MW) 
but only weakly detected within the Large Magellanic Cloud 
(LMC) and is practically absent within the Small Magel- 
lanic Cloud (SMC). The origin of the feature is generally 
believed to be aromatic carbonaceous materials, i.e. a mix- 
ture of Polycyclic Aromatic Hydrocarbon (PAH) molecules, 
which are abundant in the MW (Draine 2003; Wang et al. 
2005). A few reports of the detection of a strong 2175 A fea- 
ture in individual absorbers have been made, e.g. Srianand 
et al. (2008) and Jiang et al. (2010) analysing two individ- 
ual Mgll systems each, Zhou et al. (2010) analysing a sin- 
gle extremely dusty object, and very recently Jiang et al. 
(2011), who report 39 detections. On the other hand, statis- 
tical studies of both Call (Wild & Hewett 2005a) and Mgll 
absorbers (Wang et al. 2005; York et al. 2006) favour SMC- 
or LMC-like extinction curves. Such studies to date have 
been characterised by modest sample sizes and the E(B — V) 
values of the absorbers are typically small, which means the 
limits on the presence of any sub-population of absorbers 
possessing MW-like extinction curves are not strong. 

In this paper we employ a new, large, sample of Mgll 
absorbers, derived using quasars contained in the SDSS 
DR6-release, to investigate the form of the rest-frame ex- 
tinction curve associated with the absorber population with 
unprecedented accuracy. Adopting the approach of Wild et 
al. (2006), a full analysis of the presence of dust in the ab- 
sorbers allows a quantitative determination of the selection 
biases that effect the observed Mg II absorber sample. With 



such a quantification in hand we then review the implica- 
tions of the selection biases on i) the redshift evolution of 
the dust content of Mgll absorbers, ii) the redshift distri- 
bution of very strong Mg II absorbers, and iii) the predicted 
number of Mgll absorbers along GRB-sightlines. 

The outline of this paper is as follows. In Section 2 we 
describe the quasar and Mgll absorber catalogues which 
are used in Section 3 to estimate statistically the amount of 
dust in individual Mg II systems. We also assess the effect of 
a dust bias on the completeness of the sample in Section 4. 
In Section 5 we show results of an investigation to identify 
the 2175 A bump in the strongest Mgll systems and we look 
at the correlations between absorber EW and redshift. We 
also discuss the implications of the results for Mg II absorber 
statistics and the discrepancy in the observed frequency of 
absorbers along GRB- and quasar-sightlines in Section 5, be- 
fore summarising our conclusions in Section 6. Throughout 
this work we assume a cosmology of — 0.3, Q,a — 0.7, 
and h 0.7. 



2 QUASAR AND Mgll ABSORBER SAMPLES 

The base quasar sample consists of 91 665 objects, includ- 
ing 77 392 quasars in the Schneider et al. (2007) DR5 cata- 
logue that are retained in the later DR7 quasar catalogue of 
Schneider et al. (2010). A further 13 081 objects are quasars, 
present in the additional DR6 spectroscopic plates, identi- 
fied by one of us (PCH) using a similar prescription to that 
employed by Schneider et al. (2007), all of which are present 
in the Schneider et al. (2010) catalogue. An additional 1192 
objects, which do not satisfy one, or both, of the emission 
line velocity width or absolute magnitude criterion imposed 
by Schneider et al. (2007), are also included. None of the 
results in the paper depend on the exact definition of the 
'quasar'-sample used. 

The spectra are all processed through the sky-residual 
subtraction scheme of Wild & Hewett (2005b), resulting in 
significantly improved S/N at wavelengths A>7200A and 
consequent increased detection efficiency for Mg II absorbers 
at z>1.6. 

The quasar sample is restricted to the redshift interval 
0.5^2^3.5 and to objects with Galactic extinction corrected 
SDSS i-band magnitudes rrii ^ 19.1, leaving 52 112 quasars. 
Broad Absorption Line (BAL) quasars identified by Gibson 
et al. (2009) or from our own BAL catalogue (Allen et al. 
2011) are also excluded, leaving 48 587 quasars. 

The presence of extended wavelength intervals without 
data in the SDSS quasar spectra is not ideal for the ex- 
tinction estimation (Section 3). The sample is therefore re- 
stricted to essentially 'complete' spectra by the requirement 
that the number of valid pixels, NG00D>3500, eliminating a 
further 490 quasars. 

The final sample of 48 097 quasars is searched for Mg II 
AA2796.35, 2803.53 absorption doublets in the redshift inter- 
val 0.4^z a bs^2.2. A 'continuum' is defined for each quasar 
via the application of a simple 61-pixel median filter. The 
'difference' spectrum, to be searched for absorption fea- 
tures, is then obtained by subtracting the continuum from 
the original quasar spectrum. The absorption line search 
uses a matched-filter technique (e.g. Hewett et al. 1985) 
with two template Gaussian doublets of the appropriate 
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Figure 1. Left-hand panel: the equivalent width (EW of the Mgll A2796 line) distribution of absorbers in a sample of 19 315 systems 
from SDSS. The vertical dashed line at 1 A corresponds to lower EW-limit of absorbers retained in our sample (9719 systems). The 
sample is highly complete down to EW= 1 A. Right-hand panel: the redshift distribution for systems with EW^ 1 A. 



wavelength separation and full width at half maximum 
(FWHM) = 160 km s -1 (the resolution of the SDSS spec- 
tra), 200 km s -1 and 240kms _1 . Three values of the dou- 
blet ratio are incorporated into the search: 2:1 (correspond- 
ing to unsaturated lines on the linear part of the curve of 
growth), 1:1 (for saturated lines on the flat part of the curve 
of growth) and an intermediate case, 4:3. A further template 
with doublet ratio 1:1 and flat-bottomed absorber profiles, 
380kms _1 in extent, is also used to optimise the detection 
of very high EW doublets. At each pixel the template giv- 
ing the minimum %2 value is determined and candidate ab- 
sorbers selected by applying a threshold value of S/N>7cr. 

A total of 19 315 Mgll systems are indentified but we 
confine our analysis to the sub-sample of 9 719 systems 
with EW of the 2796 A line (W^ 2796 )^lA. The left-hand 
panel of Fig. 1 shows the distribution of W A2796 with the 
V^ A2796 =lA-limit indicated. The right-hand panel of Fig. 1 
shows the redshift distribution for the ^1 A sample. Provid- 
ing that any dust associated with the Mgll absorber pop- 
ulation results in an extinction curve whose shape is not 
strongly dependent on W A2796 , the investigation presented 
here is not sensitive to the completeness of the absorber cat- 
alogue. However, the statistical properties of the absorber 
catalogue are in excellent agreement with previous work 
(Nestor et al. 2005). The catalogue is highly complete down 
to W A2796 = lA and, from a combination of visual inspec- 
tion and investigation of the properties of other absorption 
species in 'stacked' Mgll systems, contamination by false 
detections is found to be < 5 per cent. 

There are 7318 quasars with a single intervening Mgll 
absorber, 1 019 quasars with two absorbers and 116 quasars 
with more than two absorbers. Due to the increased uncer- 
tainty in the determination of the dust content of individual 
absorption systems in quasars containing multiple absorbers 
(i.e. > 2 systems), we chose to retain only quasars with single 
or double intervening Mgll absorbers. In the case of a dou- 
ble absorber we consider only the stronger system. The vast 
majority of double absorbers involve one or more low EW 
systems, however, we remove from the sample the 8/1019 
spectra that possess two absorbers with an EW> 2.5 A to 
leave a total of 8 329 systems. The results and conclusions 
derived in the paper are not sensitive to whether we use only 
single systems or single plus double systems. 



3 MEASURING THE DUST CONTENT OF 
THE ABSORBERS 

The presence of dust in an Mgll absorber affects the ob- 
served spectrum of the background quasar. The degree of 
reddening of a spectral energy distribution (SED) is de- 
scribed by an associated extinction curve, and the dust con- 
tent is parameterised by a corresponding E(B — V) 1 (e.g. 
Kinney et al. (1994)). The large and homogenous SDSS DR6 
spectroscopic catalogue provides an opportunity to obtain 
estimates of the dust content for many thousands of ab- 
sorbers. However, one must be sure that intrinsic quasar-to- 
quasar SED variations do not produce significant systematic 
biases in, or add significant scatter to, the E(B — V) deter- 
minations of the absorbers. 

The population of Mgll absorbers is reported to have 
a mean dust content of only E(B — V) ~ 0.02 — 0.03 mag 
(Wild et al. 2006; York et al. 2006), which implies only small 
alterations to the shape of the SEDs of background quasars. 
Therefore, care must be taken to construct a quasar 'con- 
trol' spectrum, which represents accurately the spectrum of 
an unabsorbed quasar. Such a spectrum then enables us to 
isolate the effect of dust in an absorber on the background 
quasar spectrum. 

3.1 The control quasar spectrum 

We have adopted an approach similar to Wild & Hewett 
(2005a) and York et al. (2006), which involves using the 
two-dimensional redshift, z, versus magnitude, m*, plane to 
identify a sample of quasars from which to construct a con- 
trol spectrum. The use of 'neighbour' quasars with similar 
redshifts in the control sample ensures that the wavelength 
coverage of control quasars is nearly identical to that of the 
target quasar, while also ensuring that the control sample 
quasars do not have systematically different (redshift de- 
pendent) SEDs. At fixed redshift the quasar apparent mag- 
nitude distribution is equivalent to the absolute magnitude 
distribution, Mi. The large density of objects in the rrii ver- 
sus z plane allows the control sample to be defined using 

1 E(B — V) values refer to the absorber rest-frame unless specif- 
ically indicated otherwise. 
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Figure 2. Left-hand panel: flux ratio as a function of wavelength for quasar SDSS J093508.37-|-271648.6, £qso — 0.94, with an Mgll 
absorber at z a b s = 0.66. A 'control' quasar spectrum, representing the spectrum of the background quasar, is overplotted in red. Right- 
hand panel: flux ratio spectrum for the z a ^ s = 0.66 Mgll absorber in the rest-frame of the absorber, obtained by dividing the quasar 
spectrum by the control spectrum in the left-hand panel. The red line shows the best fitting SMC extinction curve, which corresponds 
to an E(B - V) = 0.15. 



magnitudes close to that of the target quasar, minimising the 
effect of any systematic magnitude-dependent SED changes. 

The control spectrum was generated for each target 
quasar, by taking the median value (at each wavelength) 
of the 50 immediately brighter and 50 immediately fainter 
quasar spectra within a fixed redshift interval of Az = 0.1 
centred on the redshift of the target quasar 2 . The median 
was chosen instead of the mean, as it was less sensitive to 
flux outliers arising from intrinsic variation among quasar 
SEDs. A control spectrum is created for each of the 9 719 
quasars with associated Mgll absorbers, but only quasars 
without a detected Mg II absorber are allowed to contribute 
to the control spectra. 

Using a fixed number of quasars to define the control 
spectrum, rather than a fixed magnitude interval, means 
that the control sample is not biased due to the systematic 
mismatch between the luminosity of the target quasar and 
median luminosity of the control quasars. 



F (A) = F 10 



2.5 



(1) 



to the flux ratio spectrum, where A (A) is the extinction at 
a particular wavelength given by 



A (A) = E (B - V) 



E x - 



(2) 



where E\-v /Eb-v is given by Pei, and Rv is the ratio of to- 
tal to selective extinction. We have adopted R v — 3.0, which 
is representative of the value in the MW and SMC 3 . Fo and 
the reddening parameter E(B — V) are left as free param- 
eters in the fit, which is carried out using the non- linear 
'Levenburg-Marquardt' algorithm (Marquardt 1963). The 
wavelength intervals corresponding to the strong telluric sky 
emission lines at AA5578.5, 6301.7, and to the prominent ab- 
sorption lines 4 of Fell, AlII, Al III, Mgl, Call, which may 
be present in the absorber spectrum, are excluded from the 
fit. 



3.2 Estimation of the E(B — V) values for the 
absorbers 

The flux ratio for an absorber is calculated by dividing 
the target quasar spectrum by the control spectrum in the 
quasar rest-frame, removing, statistically, the signature of 
the quasar SED. The resulting flux ratio spectrum is then 
moved to the absorber rest-frame (Fig. 2). 

We have used a parameterisation of the MW, LMC and 
SMC extinction curves, due to Pei (1992), and fit a curve of 
the form: 



2 The 100 quasar spectra were shifted to the quasar rest-frame, 
and normalised over the wavelength interval common to all spec- 
tra. The median value at each rest-frame wavelength was then 
calculated. 



3.3 Uncertainty in the E(B — V) determinations 

In order to test the accuracy of the method for calculating 
the E{B — V) for a particular absorber 5 , we obtained esti- 
mates for the amounts of dust in a subsample of 4 550 quasar 

3 The commonly adopted values of Ry in the SMC and MW are 
3.1 and 3.2 respectively (Fitzpatrick & Massa 2007). The shape 
of the extinction curve does not depend on the value of R v . 

4 The wavelengths of masked absorption lines are: Fe 11(1527, 
1608, 2250, 2261, 2344, 2374, 2383, 2587, 2601 A), Al 11(1671 A), 
A1III(1855, 1863 A), MgI(2853A), MgII(2796, 2803A) and 
Ca 11(3935,3970 A). The mask consists of 10 A wide intervals cen- 
tred on the absorption lines. 

5 A test of the method involved calculating E(B — V) esti- 
mates for the sample of Ca Ilabsorbers studied by Wild & Hewett 
(2005a) and no systematic differences with either redshift or 
E(B — V) were found. 
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Figure 3. The E(B — V) distribution (derived by fitting each ab- 
sorber's flux ratio with an SMC extinction law) for 4 550 quasar 
spectra, with no identified absorbers, from our quasar sample 
(Section 2). The standard deviation of &e(B-V) =0.025 quanti- 
fies the uncertainty in an individual estimate of E(B — V) due to 
quasar-to-quasar SED variations (Section 3.2). 



spectra with no identified absorbers from our quasar sam- 
ple (Fig. 3). For the non-absorber quasar sample we create 
a simulated set of absorbers with z a bs and z a bs/zqso distri- 
butions which are consistent with the corresponding red- 
shift distributions in the absorber sample. We then obtain 
E(B — V) estimates for this sample using the prescription in 
Section 3.2. The distribution is centred on E(B — V) = 0.0 
indicating no detectable systematic bias in the determina- 
tion of E(B — V), with an encouragingly small uncertainty, 
cte(b-v) ~ 0.025, in the E(B — V) estimates for individual 
quasars 6 . The spread arises from intrinsic quasar- to- quasar 
SED variations. 

Notwithstanding the use of an individual control spec- 
trum for each quasar, quasars with intrinsically unusual 
SEDs lead to flux ratio spectra still dominated by the varia- 
tion among quasar SEDs. The presence of such variations in 
turn leads to poor determinations of E(B — V). We have esti- 
mated the quality of fit of each absorber's flux ratio with an 
extinction curve by using a Kolmogorov- Smirnoff- like statis- 



tic. Our test uses the cumulative difference between the nor- 
malised best-fit extinction curve and normalised absorber's 
flux ratio, to identify a systematic shape difference and thus 
poor- fit. The maximum cumulative difference, or D ma x, is 
obtained for each absorber's SMC, LMC and MW extinc- 
tion curve fit. The D ma x distributions appear to show no 
evidence of significant bias as a function of absorber red- 
shift, E(B - V), EW and quasar magnitude (Fig. 4). 

Visual inspection of quasar spectra with large values 
of D max reveals a mixture of quasars showing i) strong 
Fell emission, ii) very blue continua or hi) red, but 'flat', 
spectra without the curvature expected due to conventional 
extinction curves. The number of such quasars is small, con- 
sistent with the tails of the distribution evident in Fig. 3, due 
to the small fraction of quasars with unusual intrinsic SEDs 
(unrelated to any signature due to intervening absorbers). 

The 2cr uncertainty in E(B — V) resulting from quasar- 
to-quasar SED variations corresponds to an E(B — V) dif- 
ference of 0.050 mag (Fig. 3). Such a difference produces a 
value of D ma x=0.04 and we therefore define spectrum flux 
ratio fits with D max >0.04 as possessing a poor fit to the SMC 
extinction curve. The number of such spectra (17/8 329) rep- 
resent an extremely small fraction of the total (0.2 per cent). 
We exclude such objects from further analysis until we con- 
sider the detailed shape of the Mgll absorber extinction 
curves in Section 5.3. 

The effectiveness of the Dm ax statistic is verified by not- 
ing that the flux ratio spectra with multiple (> 2) absorption 
systems exhibit systematically larger values of D max than 
spectra with single or double absorbers. The differences are 
due to the presence of multiple extinction signatures at dif- 
ferent redshifts, providing justification for the use of single 
and double absorber lines of sight only (Section 2). 



4 DUST OBSCURATION BIAS 

Extinction of quasar light due to the presence of dust 
in intervening absorbers, leads to objects being lost from 
magnitude-limited quasar surveys (Fall & Pei 1993; Wild 
et al. 2006) . To illustrate the effect of extinction by interven- 
ing dust on our sample, the left-hand panel of Fig. 5 shows 
the observed extinction at 7500 A (the effective wavelength 
of the SDSS i-band). Additionally, extinction due to an in- 
tervening absorber leads to a degredation of the spectrum 
S/N and hence additional absorbers are lost from the sam- 
ple, as their S/N falls below the 7a Mgll absorber detection 
threshold (Section 2). 



It is necessary to modify how the quasar control sample is de- 
fined when we approach the extremes of the mi plane where there 
are not enough quasars to make up the control-spectrum using 
the scheme described in Section 3.1. In such a case we use the 
100 quasars closest to the faintest/brightest quasar at a given 
redshift. The effectiveness of the modified sample definition is 
demonstrated by calculating E(B — V) estimates for a sample of 
the 51st brightest quasars using control spectra defined using both 
the ±50 neighbouring quasars and the next 100 fainter quasars. 
Average differences in the E(B — V) estimates of only 0.001 at 
the faint end and 0.0025 at the bright end, demonstrate that any 
systematic bias in the E(B — V) estimates is insignificant. 



4.1 Extinction effect 

We have used the magnitude and redshift distribution of the 
quasars in the quasar sample (Section 2) to obtain a prob- 
ability that a quasar behind a dusty absorber will be lost 
from the sample. The presence of an absorber in a quasar 
spectrum leads to an extinction of quasar light Ai (Fig. 5, 
left-hand panel). If the quasars in the sample 7 were placed 
behind the absorber, objects fainter than 19.1 — Ai will be 



7 We consider only absorbers that are displaced at least Az>0.1 
from the target quasar 
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Figure 5. Left-hand panel: the observed (i-band) extinction caused by intervening LMC-type dust as a function of absorber redshift. 
The different curves are for different values of E(B — V). Right-hand panel: the i-band magnitude versus quasar redshift for all objects 
in the quasar sample (Section 2). Objects to the right of the vertical red line are quasars which are Az>0.1 more distant than an 
absorber (at z a t, s = 1.0), and the horizontal dashed line corresponds to the detection limit of quasars placed behind the absorber with 
E(B — V) = 0.125. The extinction, (A{), is calculated as illustrated in the left-hand panel. 



lost below the vra — 19.1 flux limit. We can therefore quan- 
tify the probability an absorber of specified E(B — V) and 
redshift is lost from the sample via 



iVqso (< 19.1 - Aj, z ahs < z qso - 0.1) 
iVqso « 19.1, Z ahs < Zqso ~ 0.1) 



(3) 



where N qso is the number of quasars satisfying the condition 
in brackets. An illustration of the bias calculation is shown 
in the right-hand panel of Fig. 5 for an absorber at z a bs = 1-0 
with E(B — V) = 0.125. 



sorber will be lost from the sample due to the reduction in 
the spectrum S/N, 



PS/N = 1 



(S/N 



S/N q 



S/N q 



^ 7a) 



(S/N a 



(5) 



where S/N qgo is the S/N of a given quasar spectrum in the 
base quasar sample. Ps/n is calulated using only the quasars 
that remain in the sample following the determination of 
Pdob in Section 4.1. 



4.2 Signal-to-noise ratio effect 

For absorbers that remain within the sample there is a fur- 
ther reduction in the probability of detection due to the 
decrease in the S/N of the quasar spectrum as the quasar 
becomes fainter because of the extinction from dust. Taking 
only the quasars that remain above the m— 19.1- Ai thresh- 
old determined above, we can determine how a quasar spec- 
trum's S/N degrades for a given extinction Ai by examining 
the median S/N versus rrii for the quasar spectra. Consider 
the effect of placing an absorber of specified E(B — V) and 
redshift in front of a quasar with a given ra*, and hence S/N. 
The extinction due to the absorber makes the quasar fainter 
and the Mgll absorber may no longer be detectable given 
the reduced S/N of the spectrum. An Mgll absorber will be 
lost from the sample when 



S/N a 



^/^qsoext n 

x ~ — rr"~ ^ 7a 



S/N, 



(4) 



qsoabs 



holds true, where S/N abs is the S/N of the Mgll absorption 
feature, S/N qsoext is the degraded S/N of the quasar spec- 
trum, S/N qsoabs is the S/N of the absorber's target quasar 
spectrum, and 7a is the Mgll absorber detection threshold. 
We can therefore obtain a separate probability that an ab- 



4.3 Correcting the statistics 

The base absorber sample is highly incomplete at large 
values of E(B — V), and we provide a correction to the 
E{B — ^-distributions for both the extinction and S/N ef- 
fects described in this section. Each absorber is assigned an 
E(B — V)- weighting due to both obscuration effects accord- 
ing to 

W = 1-Pdob (E(B-V),z) X 1-Ps/n (E(B-V),zY (6) 

which allows the observed E(B — ^-distributions to be 
corrected. The weighting, w, becomes undesirably large for 
absorbers with large E(B — V), where the probability an 
absorber is included in the sample is small. We therefore 
impose a conservative upper limit to the value of u>=6.67, 
equivalent to the probability an absorber is included in the 
sample of 0.15. The multiplicative prescription for calculat- 
ing w is valid as we only apply the S/N corrections to the 
population of absorbers which are not lost due to dust ex- 
tinction effects. Thus, suppose we have a sample of absorbers 
from which 30 per cent of quasars fall below the magnitude 
limit due to extinction. The reduction in the S/N is calcu- 
lated for the remaining 70 per cent of quasars and (in the 
case that the S/N effect reduces the number of absorbers 
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Figure 4. The Kolmogorov-Smirnoff-type D ma x statistic for the 
base sample of Mg II absorbers, fitted with an SMC-type extinc- 
tion law, as a function of (top to bottom panels) absorber redshift, 
E(B — V), equivalent width and quasar magnitude. The median 
points show no significant bias in the quality of fit as a func- 
tion of any of the parameters. The horizontal line indicates the 
D m ax=0.04 threshold, used to define systems with a 'poor fit' to 
the SMC extinction law. 

by 10 per cent) the fractional completeness of the absorber 
sample is 0.7 x 0.9 = 0.63. 

4-3.1 Accuracy of the corrected E(B — V) estimates 

The observed E(B — V) distributions show a spread at fixed 
EW due in part to quasar- to- quasar SED-variations, quan- 
tified in Section 3.3 and illustrated in Fig. 3. In principle, 
a full deconvolution of the SED-induced variations could be 
undertaken using the observed E(B — V) distributions to 
recover the intrinsic E(B — V) distributions. However, such 
a procedure is both involved and unstable in portions of the 
E(B — V) versus EW plane. 

The use, in Section 5, of the median E(B — V) values 
for the absorber population as a function of absorber EW 
already mitigates the impact of the intrinsic quasar SED- 
induced variations. However, the procedure described in the 
previous sub-sections does result in a small overestimation 
of the population median E(B — V). The worst-case situa- 
tion occurs when considering a population of absorbers with 
no dust, subject to the quasar SED-induced spread. A his- 
togram of 'observed E(B — V) \ exactly as shown in Fig. 3 is 
thus obtained. The observed median E(B — V) is only a few 
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Figure 6. Left-hand panel: Dust obscuration bias for SDSS 
quasar absorption systems as a function of absorber redshift 
( z abs) an d E(B — V). The contours represent the probability that 
an individual absorber will not be observed. The dashed horizon- 
tal line represents the E(B — V) = 0.15 value which corresponds 
to a dust obscuration bias of 0.85 at z a b s = 2. In Section 5.2 we 
use this E(B — V) value to define a conservative E(B — V)-cut 
to examine any statistical trends with redshift 



thousandths of a magnitude from zero but, following cor- 
rection of the E(B — V) distribution according to Equation 
6, the median E(B — V) becomes +0.016 mag. The origin 
of the positive value derives in part to the small asymme- 
try to positive values in the histogram but (mostly) to the 
application of the correction scheme only to objects with 
apparently positive E(B — V) values. 

The overestimation of the median E(B — V) is small 
under worst-case conditions and becomes much smaller 8 as 
populations with dust-induced E(B — V) are considered and 
the fraction of absorbers scattered artificially to negative 
E(B — V) drops. Our conclusions in the paper, which derive 
from the presence of median E(B — V) values > 0.07 mag 
for absorbers with high Mg HEW, are essentially unaffected 
by the form of correction adopted to account for the small 
overestimation of the median E(B — V) values. Results ob- 
tained after applying a uniform correction of -0.016 mag to 
all the median estimates 9 and also no correction at all are 
not significantly different. We thus adopt a simple scheme 
that closely approximates the actual situation, applying a 
correction of -0.016 mag to the median E(B—V) estimates of 
the absorber samples with MgllEW ^ 3.0 A and no correc- 
tion to the median E(B-V) estimates for Mg HEW > 3.0 A. 



8 The offset is < +0.005 for the median E(B-V) estimate of the 
absorber sample with MgllEW > 3 A. 

9 This correction is performed by subtracting 0.016 from each 
E(B — V) value in the absorber sample and recalculating the 
weighted median. 
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Figure 7. E(B — V)smc as a function of Wq 2796 for a sample 
of Mgll absorption systems with accurate E(B — V)smc values 
(i.e. those with D^^. < 0.04). The blue points represent the me- 
dian values for a constant equivalent width bin of 0.7 A, and the 
red points represent the median values adjusted for dust obscura- 
tion bias. The error bars are calculated by using the 'bootstrap' 
method, and the solid red line represents the best fitting power- 
law to the bias-corrected data. At low values of Wq 796 there is 
little evidence for significant E(B — V) associated with the ab- 
sorbers and the distribution of E(B — V) values is dominated by 
the quasar SED-induced spread illustrated in Fig. 3. The true un- 
certainty in the estimation of the median E(B — V) at low values 
of Wq 796 is limited by systematics at the level of ~ 0.005. 



5 RESULTS 

The amount of dust in Mgll absorption systems is gener- 
ally found to increase as a function of the Mg II EW (York 
et al. 2006; Menard et al. 2008). This increase has impor- 
tant implications for the understanding of galaxy evolution, 
as EW for the saturated Mg II absorbers measures the veloc- 
ity spread in the gas (York et al. 2006) and there is a strong 
dependence of EW on the level of associated star- format ion 
activity (Menard et al. 2009). 



5.1 Dust and absorber Equivalent width 

The accuracy of the E(B — V)smc 10 estimate for each ab- 
sorber was ensured by removing all spectra with a poor fit 
to the SMC extinction curve, i.e. > 0.04 (Section 3.3). 

We chose the SMC curve as the average reddening proper- 
ties of Mg II absorbers are found to be well characterised by 
an SMC-type extinction law (see Section 5.3). 



10 To make the form of the extinction curve used to obtain 
E(B — V) values explicit, where relevant, we use a subscript 
'LMC, 'MW, 'SMC',... in the remainder of the paper 



5.1.1 E{B - V) Distributions 

The dependence of dust content on Mgll EW is shown in 
Fig. 7. The median E(B — V)smc is calculated for absorbers 
with EWs in the intervals 1.0-1.7 A, 1.7-2.4 A, 2.4-3.1 A, 3.1- 
3.8 A, 3.8-4.5 A, and 4.5-5.0 A. The equivalent median E(B — 
V^smc values incorporating the bias corrections in Equation 
7 are also shown. The error bars have been calculated using 
the 'bootstrap' method but, at low W^ 796 , the errors are 
limited by systematic uncertainties at the level of ~ 0.005. 
The line in Fig. 7 corresponds to the best fitting power-law 
to the median bias-corrected E(B — V)smc values of the 
form 

(E(B - V)smc) = A (W X2796 ) a , (7) 

where A = (8.0 ± 3.0) x 10" 4 and a = (3.48 ± 0.3) and the 
model is applicable for 1.0<Wo^ 5.0 A. 

The number and distribution of objects evident in Fig. 7 
with E(B — V)smc < are consistent with the expected 
observational uncertainty in the E(B — V) determinations 
(Section 3.3). A series of E(B — V)smc distributions for dif- 
ferent bins of mean Mgll EW 1.2 A, 2.0 A, 2.8 A and 3.6 A, 
are shown in Fig. 8. 

5.1.2 Comparison with other work 

Menard et al. (2008) provide a parameterisation of the 
E(B — V) content of strong Mgll absorbers as a function 
of EW that has started to see wide use in the literature 
(e.g. Elfasdottir et al. 2009; Sudilovsky et al. 2009). Menard 
et al.'s investigation presented observed-frame E(B — V)smc 
distributions, which contrasts with the absorber rest-frame 
parametrizations, made possible by the direct use of the 
SDSS quasar spectra, presented here. 

The observed- frame E(B — V)smc dependence on Mgll 
absorber EW (equation 7 of Menard et al.) can be related 
to our rest-frame E(B — V)smc dependence via 

E(B - V) ohs = E(B - V%st (1 + z) 1 ' 2 . (8) 

Fig. 9 shows our data in the E(B — V) bs versus EW 
plane (c.f. fig. 9 of Menard et al.)) along with the median 
E(B — V)obs points, median dust-corrected E(B — V)smc 
values and the power- law fit (Equation 7), all transformed 
into the observed frame. There is good agreement between 
the Menard et al. observed median E(B — V)smc values and 
those determined here. However, our median dust-corrected 
E(B — V)smc values lie significantly above Menard et al.'s 
observed values. The large differences can be ascribed to 
the ability to quantify the effects of E(B — y)-bias on the 
observed Mgll sample, set out in Section 4. These dust ob- 
scuration effects while linked to EW, also become more pro- 
nounced at higher redshift. 

5.2 Redshift evolution of E(B - V) 

The ability to constrain how the dust content of Mgll ab- 
sorbers varies with cosmic time, from 3.3 Gyr at z — 2.0 
to 8.6 Gyr at z — 0.5, when a dramatic reduction in in the 
overall star formation rate density of the Universe (Hopkins 
2004; Menard et al. 2009) occured, provides a significant 
constraint on models of galaxy evolution. 
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Figure 8. The E(B — V)smc distributions of Mgll absorption systems for different equivalent width ranges. The solid line represents 
the distribution of observed E(B — V)sMC- es timates, and the dotted line corresponds to the distribution corrected for bias (Section 4). 
The dotted lines represent estimates of the intrinsic E(B — V)smc distribution, which exhibit more pronounced high-E(B — V) tails. 



Notwithstanding the quantitative corrections to the ob- 
served E(B— V)smc distribution of the Mg II absorbers (Sec- 
tions 4 and 5.1) the sample suffers from significant incom- 
pleteness for absorbers with large E(B— V) at high redshifts. 
To investigate the evolution of the absorbers as a function 
of redshift we therefore confine our analysis to systems with 
E(B - V)smc < 0.15. The adopted limit to E(B - V)smc 
corresponds to that of an absorber Pdob>0.85 n at the max- 
imum absorber redshift, 2=2.0, employed (see Fig. 6). 

5.2.1 Results and comparison with previous work 

The restricted absorber sample (E(B — V)smc < 0.15) was 
divided into three redshift bins, with average redshifts of 
0.77, 1.25, and 1.68. The resulting bias-corrected median 
rest-frame E(B — V)smc are fit with a power- law of the 
form 

E(B - V) oc (1 + z) a , (9) 

where a — — 0.2 ± 0.3, which shows no significant evolution 
with redshift. This determination is consistent with Menard 

11 Application of the bias-correction to absorbers with Pdob ^ 
0.85 can cause undesirable divergence in the population statistics. 



et al. (2008) who find a modest evolution of dust content 
with redshift (a = -1.1 ± 0.4). 

While this analysis of the bias-corrected median points 
has deliberately been confined to absorbers with E(B — 
V) smc < 0.15 it is also possible to test whether the entire ob- 
served absorber sample is consistent with no evolution as a 
function of z a bs- Taking the bias corrected E(B — V)smc dis- 
tribution of absorbers in the low-redshift (z a bs) =0.77 slice 
we calculate the observed E(B — V)smc distribution, by ap- 
plying the inverse of the bias corrections to the (z a bs) =0.77 
intrinsic distribution, centred on the higher redshift slices 
(^abs) = 1-25 and (z a bs) = 1.68 respectively (Fig. 10). The 
resulting distributions for the higher redshift slices are con- 
sistent with the observed histograms, confirming the lack 
of evidence for significant evolution in the distribution of 
E(B — V)smc with redshift. 



5.3 Nature of the dust 

Determining the nature of dust in Mg II absorption line sys- 
tems is important as it provides a way to constrain the 
chemical evolution of galaxies over a range of cosmic time. 
Of particular interest is whether the strong ~ 2175 A fea- 
ture observed in the spectrum of the MW is present in 
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Figure 9. The observed- frame E(B — V)smc -distribution of 
Mgll absorbers as a function of Wq 2796 . The observed- frame 
E(B — V)smc values (calculated using Equation 8) are plot- 
ted to enable comparison with the Menard et al. (2008) result 
shown as the solid-green line. The blue points are the median 
observed- frame E(B — V)smc values calculated using EW bins 
of AEW=0.7A, and the red points represent the median values 
adjusted for dust-obscuration bias. The error bars are calculated 
by using the 'bootstrap' method, and the red line represents the 
best fitting power-law to the bias-corrected data. The true uncer- 
tainty in the estimation of the median E(B-V) at low values of 
Wq 796 is limited by systematics at the level of ~ 0.005. 

our Mgll absorber sample. In general, the SMC extinction 
curve, rather than that of the MW, is found to best de- 
scribe the average reddening properties of Mgll absorbers 
(Menard et al. 2005; York et al. 2006; Menard et al. 2008) 
with low E(B — V). Taking advantage of the large sample of 
absorbers, and specifically the availability of a number of ab- 
sorbers with significant E(B — V), we investigated the form 
of the extinction curve via construction of a composite spec- 
trum and consideration of the properties of the absorbers 
with the very highest estimates of E(B — V). 

5.3.1 An optimal high-E(B — V) composite absorber 
spectrum 

Despite claims of being able to see evidence for the 
2175 A bump in single systems (Srianand et al. 2008), in- 
dividual extinction curves have a low S/N. Using the much 
larger sample of absorbers presented here, we can combine 
a large number of extinction curves statistically, using an 
approach similar to Wild & Hewett (2005a), to improve the 
S/N and increase the sensitivity to the presence of any fea- 
tures. 

An important consideration is the choice of which sys- 
tems to combine. The optimum number of spectra to co-add 
is a trade-off between the number of spectra (in principle the 
more the better), and the median E(B — V) of the coadded 



Figure 10. The normalised distribution of absorbers for the full 
sample (in different redshift bins) as a function of E(B — V)smc- 
The histograms have been normalised to enable comparison of 
redshift bins containing different numbers of absorbers. The black 
solid line corresponds to the distribution of absorbers in the low- 
redshift ((z) = 0.77) bin when corrected for dust bias effects. 
The green and red histogram are the observed distributions (i.e. 
uncorrected for dust bias) of E(B — V) for the higher redshift 
bins ((z) = 1.25 and (z) = 1.25) respectively. The dotted lines 
represent the result of applying the dust obscuration effects to 
the black histogram at these higher redshifts. The consistency 
between the observed histograms and the predicted distributions 
again highlight that there is no evidence for significant evolution 
of the dust content of Mgll absorbers as a function of redshift. 



spectra (which favours a smaller number of spectra with the 
largest E(B - V». 

We also need to consider the goodness-of-fit of the spec- 
tra that will make up the composite spectrum, and so it was 
necessary to make a cut based upon the D m ax statistic (Sec- 
tion 3). However, it is conceivable that a system with a large 
SMC-D m ax value, and thus poor SMC fit, will in fact be 
well characterised by a MW extinction curve. We certainly 
do not want to remove such objects which exhibit a MW 
type dust signature from the composite sample. This D max 
discrepancy (between SMC and MW extinction curves) be- 
comes more significant at higher E(B — V) values, and we 
can quantify, for a given absorber redshift, the E(B — V)smc 
value which corresponds to a poor- fit D ma x = 0.04. We have 
determined how this E(B — V)smc cutoff varies as a function 
of redshift and the results are shown in Fig 12. 

Absorption systems in the shaded region of Fig 12, are 
those which satisfy the parameterised condition: 

E(B -V)smc >ae~ bz + c (10) 

where a = 38.7, b = 4.55, and c = 0.04. For such absorbers, 
the form of the extinction curve, SMC or MW, can be de- 
termined via their SMC/MW D max statistics. 
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Figure 13. Top panel: Wavelength versus flux ratio for a composite mean spectrum created from 211 absorber flux ratios which 
satisfy Equation 10. Additional composites of each absorber's best fit SMC, LMC and MW extinction curves are also shown with mean 
E(B — V)smc — 0.12, E(B — V)lmc — 0-16 and E(B — V)mw — 0-13. The SMC curve is the best fit to the composite flux ratio spectrum. 
The black box defines the zoomed-in section around the 2175 A bump. Bottom-left panel: The residuals between a mock coadded curve 
(containing x per cent LMC spectra, y per cent MW spectra and the rest SMC spectra), and the composite flux ratio. The best-fit 
corresponds to a composite containing only SMC and LMC extinction curves in the ratio 65:35. Bottom-right panel: The zoomed-in 
section of the composite spectrum with green, and red solid lines corresponding to the SMC and LMC extinction curve composites 
respectively. The gold-dashed line corresponds to the coadded curve containing SMC spectra and LMC in the ratio 65:35, which best 
characterises the composite flux ratio. 



5.3.2 Results from the composite spectrum 

We created a composite spectrum (Fig 13) from 211 ob- 
jects satisfying the condition in Equation 10, by taking the 
mean 12 flux value at each wavelength. Three further com- 
posites were created from the fitted extinction curves (SMC, 
LMC, MW) corresponding to each absorber. This allows 
for the characterisation of the type of dust by investigat- 
ing whether the stacked spectrum is well fit by the SMC, 



12 The mean (rather than the median) is now used to generate 
the composite spectrum to allow us to estimate the contribution 
of SMC-/LMC-/MW-type extinction to the composite. 



LMC or MW curves. The resulting curves possess E(B — V) 
values of 0.12, 0.16, and 0.13 for the SMC, LMC and MW 
respectively. As predicted, Fig 13 shows that on average our 
carefully-selected high E(B — V) absorber subsample ex- 
hibits an extinction curve very similar to that of the SMC. 

We can quantify the significance of the SMC- versus 
LMC- and MW-type reddening by examining the relative 
contributions of SMC-, LMC- and MW-type dust to the 
composite spectrum in the region around the 2175 A fea- 
ture (Fig. 13- {bottom-right)) . We obtain an estimate for the 
ratio of SMC-to-LMC-to-MW type spectra which contribute 
to our composite, by coadding a number of SMC, LMC, and 
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Figure 11. The E(B — V)smc of Mgll absorbers as a func- 
tion of absorber redshift. The green square and blue diamonds 
correspond to the Zhou et al. (2010) and Srianand et al. (2008) 
objects respectively. The red crosses represent the objects most 
likely to contain signatures of the 2175 A bump, which make 
up the LMC coadd sample in Fig. 15. The horizontal line at 
E(B — V)smc = 0-15 is included to show which absorbers make 
it into the redshift evolution sample described in Section 5.2, and 
the line at E(B — V)smc = 0-30 shows which objects make it into 
the high-E(B — V) sample described in Section 5.3.3. 

MW extinction curves and minimising the residuals 13 be- 
tween the coadd and the composite spectrum (in the range 
1900-2500 A). Formally, we found that a combination of 137 
SMC-, 74 LMC- and MW-curves provides the best fit to 
the composite spectrum in the 2175 A feature region. The 
best-fitting combination is shown in Fig. 13- {bottom-right) . 
The analysis shows that of the absorbers which satisfy Equa- 
tion 10, their extinction curves are well reproduced by a mix 
of SMC and LMC extinction curves in the ratio 65:35 with 
an uncertainty of ~±10 in the ratio. 

The existence of absorbers with LMC-type dust can be 
demonstrated more directly via consideration of Fig. 14, 
which shows the relationship between (x-axis) and 

both L> m ^ x c and (y-axis) for the 211 absorber sys- 

tems, in general, D^ C < Dm^, which is to be expected if 
the majority of absorbers contain SMC- or LMC-type dust. 
Coadding the flux ratio spectra for the subset of 18 absorbers 
with > 0.02 and L> m ^ x c < 0.02 produces a composite 

flux ratio (Fig. 15) the shape of which is closely reproduced 
by the LMC extinction curve. The mean E(B — V)'s for the 
SMC, LMC, and MW type extinction curves in this compos- 
ite sample are 0.16, 0.21, and 0.18 respectively. It can be ar- 
gued that the result must be true by construction, with indi- 
vidual objects chosen to possess LMC-like extinction curves. 
However, the procedure is no different from the identifica- 

13 The squared deviations between the underlying composite 
spectrum and the coadded extinction curves. 



Figure 12. The E(B — V)smc value which corresponds to a poor- 
fit Dmax = 0.04 statistic between an SMC and MW extinction 
curve, as a function of redshift. The solid black line corresponds 
to a parameterised fit to the black points, and the shaded region 
above the line represents systems which can, in principle, be dis- 
tinguished as containing either SMC or MW type dusts on the 
basis of their SMC/MW Dmax statistic. 

tion of individual systems (Srianand et al. 2008, e.g.) from 
among large samples of quasar spectra (Srianand et al. 2008, 
e.g.). More convincingly, the 18 systems represent a signifi- 
cant minority of the base sample of 211. Furthermore, only a 
tiny fraction of absorption systems were eliminated from the 
sample via a D m ax threshold (Fig. 4) and it is very unlikely 
that unrelated intrinsic quasar- to- quasar SED variations are 
combining to reproduce a 2175 A feature and overall shape 
that reproduces the form of the LMC-extinction curve to 
the accuracy shown in Fig. 15. The objects contributing to 
the composite are indicated in Fig. 11. With only half of 
these objects present in the restricted sample used to inves- 
tigate redshift evolution (Section 5.2), it is not possible to 
say anything useful about the strength of the 2175 A bump 
with redshift. 

York et al. (2006) perfomed a statistical analysis of 809 
Mgll absorption systems from SDSS DR1, and searched for 
the 2175 A feature in the dustiest 111 systems with an ob- 
served frame colour-excess 14 A (g — i) > 0.2. They find an 
average E(B — V)smc = 0.0805 and no evidence of a MW 
dust signature in their composite spectrum. They also pro- 
vide an estimate of the fraction of lines of sight that could 
have MW-type extinction, and provide an upper limit of 30 
per cent. 

Although we have reached a similar conclusion regard- 
ing the nature of dust, our sample of Mgll absorbers is an 
order of magnitude larger, and we have created a composite 

14 This excess corresponds to the difference between the actual 
colours of a quasar and the median colours of a quasar at that 
redshift. 
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Figure 14. -Dmax for the SMC curve versus the D ma x for the 
LMC curve (red-points) and the MW curve (blue-points). The 
absorbers plotted are the 211 where discrimination between dif- 
ferent extinction curves is possible according to the condition of 
Equation 10. The dashed lines represent the average poor-fit cut- 
off (Section 4) of Dmax = 0.04. The circled objects are the two 
absorbers discussed by Srianand et al. (2008), which are found 
to contain signatures of the 2175 A bump, and the solid box in- 
dicates objects which are very likely to contain bump signatures 
based upon their LMC/SMC D ma x statistic. 



spectrum with a much larger average E(B — V)smc- The 
sample used has also been carefully constructed to ensure 
that only objects with accurate estimates of E(B — V) (i.e. 
those with < 0.04) are allowed to contribute to the 

composite. We provide a tighter constraint on the fraction 
of lines of sight with (the much less extreme) LMC-type ex- 
tinction, as 35±10 per cent, compared to York et al.'s upper 
limit of 30 per cent with MW-type extinction. 

A recent study by Jiang et al. (2011) produced a cat- 
alogue of 39 candidate strong Mgll absorbers, with EW> 
lA and 1.0 < z < 1.86, showing evidence for the pres- 
ence of a weak 2175 A feature. We find 19 of the Jiang 
et al. (2011) objects in our sample. All possess moder- 
ate E(B — V)smc °° 0.1, and 15/19 satisfy the condition 
£>max < Jiang et al. (2011)'s highlighting of these 

absorbers adds further weight to the conclusion that a sig- 
nificant fraction of absorbers with intermediate E(B — V)s 
show evidence for the presence of a weak 2175 A feature. 



5.3.3 The dustiest absorbers and the shape of the 
extinction curve 

Very few absorbers at cosmological distances with un- 
ambiguous 2175 A features have been discovered. As dis- 
cussed above, Srianand et al. (2008) find two Mgll ab- 
sorbers at z^l.3, which are best-fit by LMC-type extinc- 
tion, with E(B — VJlmc's of 0.27 and 0.36 respectively. 
Our analysis also finds these absorbers to be best-fit by 



LMC-type extinction, with corresponding E(B — V)lmc's 
of 0.30 and 0.36, consistent with their result. The two ob- 
jects are present in our composite subsample and are in- 
dicated in Figs. 11 and 14. More recently, Jiang et al. 
(2010) find two Mgll absorbers at z^l.4, which may show 
detectable 2175 A features. Our analysis finds the spec- 
trum of SDSS J012147.73+002718.7 to be well fit by an 
LMC extinction curve with E(B — F)lmc = 0.17, whereas 
the other spectrum (SDSS J145907.19+002401.2, hence- 
forth J1459+0024), is best fit by the SMC curve with 
E(B - V)smc = 0.31. Jiang et al. find that J1459+0024 
is not particularly well fit by any SMC, LMC or MW mod- 
els because the extinction bump is much broader than that 
of the average LMC extinction curve. A detailed discussion 
of the shape of the extinction curve for J1459+0024, and 
the extinction curve of other high E(B — F)smc objects is 
presented below. 

Focussing on absorbers with comparable, or greater, 
estimates of E(B — V), assuming an SMC-type extinction 
curve, there are 25 systems with E(B — V)smc^0.30. Un- 
surprisingly, the majority of absorbers have relatively low 
redshifts, £ a bs<l, although a small number extend out to 
nearly 2 a bs=l-5. Two of the systems are present in the sam- 
ple of absorbers with detectable Call absorption discov- 
ered by Wild & Hewett (2005a) and Wild et al. (2006). The 
most extreme system, with E(B — V a )smc=0.60, has been 
the subject of an earlier study. Zhou et al. (2010) 's analysis 
of SDSS J100713.68+285348.4 highlights the strength of the 
2175 A feature in the strong Mgll absorber, at z a bs=0.884, 
which also shows extremely strong Call absorption, al- 
though they do not note the connection to the Ca II absorber 
population. Employing constraints on the shape of the ex- 
tinction curve extending to the rest-frame near- infrared, pro- 
vided by JHK photometry from 2MASS (Skrutskie et al. 
2006), Zhou et al. (2010) present an extinction curve us- 
ing the parametrization of Fitzpatrick & Massa (2007). The 
form of the curve includes a very strong 2175 A feature, su- 
perposed on a rather greyer background compared to any 
of the SMC, LMC or MW extinction curves. Zhou et al. 
(2010) 's fit to the SDSS spectrum and the broadband pho- 
tometry leads to a smaller value of E(B — y)=0.28 but with 
a relatively large value of itV=3.87. 

Using our full sample of 25 objects with E(B — 
V) smc ^0.30, it is possible to gain further insight into the 
shape of the extinction curve associated with the dustiest in- 
tervening absorbers. Table 1 provides a summary of informa- 
tion for the 25 systems and the background quasars, includ- 
ing Y JHK photometry from the UKIDSS survey (Lawrence 
et al. 2007) where available. The population is not signifi- 
cantly influenced by the effects of extinction from two or 
more absorbers along the line of sight to a quasar. An ex- 
tensive search for strong absorbers, including the Fell line 
complexes at A2344.2 and A2586.7, covering the redshift in- 
terval 0.36<z<2.88, shows six of the quasars possess two 
absorbers with Wo A2796 >1.0 A, whereas three are expected 
by chance. 

A striking feature of the absorption line spectra is the 
almost ubiquitous presence of Call AA3934.8,3969.6 absorp- 
tion. Fourteen of the 23 absorbers with z a bs<l-31, where 
Call lies within the SDSS spectra, individually show evi- 
dence for Call absorption and a median composite of the 
23 absorbers shows a rest-frame EW of Call A3934.8=0.38± 
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0.05 A. The absorption systems are thus closely related to 
the (observationally) rare, dusty Call absorber population 
identified by Wild & Hewett (2005a), the properties of which 
are discussed by Wild et al. (2006, 2007). The presence of 
strong absorption due to a species that is so strongly de- 
pleted onto dust as Ca II is most naturally explained by the 
presence of very large gas columns (Wild et al. 2006) and our 
interpretation is that these absorbers will possess hydrogen 
columns as large as DLAs. 

Using a model unreddened quasar spectrum (Maddox 
et al. 2008) we essentially confirm the results of Zhou et al. 
(2010) for the shape of the extinction curve for the ab- 
sorber in SDSS J100713.68+285348.4. It is possible to un- 
dertake an identical analysis for the six additional absorbers 
that possess YJHK photometry 15 . In all cases, adopt- 
ing our E(B - V) values (using any of the SMC, LMC 
or MW curves) from the SDSS spectra produces predic- 
tions that are grossly inconsistent with the SDSS (z-band) 
to YJHK colours. Specifically, compared to any of the 
SMC/LMC/MW-extinction curves, the shape of the extinc- 
tion curve probed by the SDSS spectra of the absorbers at 
redshifts z a bs— 0.75 (rest-frame wavelengths ^2100-5000 A) 
is much steeper than the curve at rest-frame wavelengths 
^4500-14 000 A probed by the broadband i to K colours. In 
other words, the extinction curve is significantly greyer at 
> 5000 A, with local extinction curves predicting a median 
i — K colour 0.72 mag redder than that observed for the six 
absorbers. 

While the five-parameter parametrization of Fitzpatrick 
& Massa (2007) would provide accurate fits to the SDSS 
spectra and the i through K photometry of the six ab- 
sorbers, we chose to adopt the established one-parameter 
family of extinction curves presented by Fitzpatrick (1999), 
in which the shape of the untraviolet to near-infrared ex- 
tinction curves depends on the value of Rv (see fig. 7 of 
(Fitzpatrick 1999)). 

We find satisfactory fits to both the flux ratio spec- 
tra derived from the SDSS quasar spectra and the iYJHK 
broadband photometry of the six objects using a Fitzpatrick 
(1999) curve with R v =2.1 (see Appendix A). The Anax- 
statistic distribution is as good as that from any of the SMC, 
LMC, MW or Zhou et al. extinction curves. The median 
i — K colours (model - observed) using the Fitzpatrick and 
Zhou curves are +0.09 mag and -0.04 respectively (cf. the 
+0.72 value using any of the SMC, LMC or MW curves). In- 
dividual iYJHK magnitudes are reproduced to the level of 
a c± 0.12 mag using either the Fitzpatrick or Zhou curves but 
with some differences reaching 0.3 mag. The small number 
of absorbers, combined with the effects of intrinsic quasar 
variability (between the epochs of the SDSS and UKIDSS 
photometry) and object to object SED variation, preclude 
the generation of more exact constraints on the form of the 
extinction curve but the need for a significantly greyer shape 
in the near-ultraviolet and optical is clear. 



15 The faint ness of the quasars, combined with the bright limit 
to the 2MASS photometry, means that choosing to analyse just 
those quasars that are detected in 2MASS (at low S/N) results 
in a significant bias, making the detected subset artificially red in 
i — K and we thus confine our attention to quasars with UKIDSS 
photometry. 



The inferred values of E(B — V) for the six systems em- 
ploying the Fitzpatrick iiV=2.1 curve are given in Table 1. 
We note that the form of the Fitzpatrick Rv -dependent ex- 
tinction curves at low values of Rv reproduces the shape of 
the extinction curve in the near ultraviolet with a 2175 A fea- 
ture of strength comparable to that in the MW, combined 
with a somewhat steeper increase in extinction with de- 
creasing wavelength, a key feature of the behaviour found 
by Fitzpatrick & Massa (2007). From the data available, 
it is not possible to conclude whether Zhou et al. (2010) 's 
modeling, employing a super-strong 2175 A feature, or the 
itV-dependent form (with an extreme value of iiV=2.1) of 
Fitzpatrick (1999), is superior but the simple one-parameter- 
dependent form, with a 2175 A feature of constant strength, 
is attractive. In the Appendix we provide both a graphi- 
cal representation of the extinction curve and quantitative 
values of E(X - V)/E(B - V) as a table. 

The difference in the form of the extinction curve (from 
SMC/LMC/MW curves) for absorbers with high E(B-V), 
including the likely presence of a strong 2175 A feature, 
raises the question of whether there is a systematic depen- 
dence of the extinction curve shape with E(B — V). While 
the detection of the 2175 A feature in individual spectra with 
small E(B — V) is not viable, the prevalence of an SMC-like 
curve for absorbers with low extinctions (E(B — U)<0.05) is 
unambiguous; any extinction curve with a 2175 A feature as 
strong as that in the MW can be ruled out with confidence 
(Fig. 16). For absorbers with intermediate E(B — U)~0.15 
we find evidence for the presence of a weak 2175 A feature in 
a significant fraction of systems (Section 5.3.2). The results 
for systems with larger E(B — U)smc require confirmation 
using a much larger sample of absorbers with near-infrared 
photometry but for the first time there is an indication of a 
strong systematic change in the form of the near-ultraviolet 
to near-infrared extinction curve as a function of E(B — U), 
with the shape of the curve in the dustiest systems differing 
from any of the standard curves in the local Universe. 

5.4 Implications for Mgll absorber statistics 

The quantitative dust-dependent bias, formulated in terms 
of the absorber E(B — V), derived in this paper has sig- 
nificant implications for the study of the Mgll absorber 
population. Fig. 17 shows the distribution of observed EW 
and redshift for our absorber sample. Also shown are the 
bias-corrected distributions, calculated using the object by 
object bias-corrections (equation 6). A significant fraction 
of absorbers are missing from the sample, 24 ± 4 per 
cent for absorbers with EW>1.0A in the redshift interval 
0.4 < z < 2.2. The missing fraction rises to 34 ± 2 per cent 
for absorbers with EW>2.0 A. The formal errors on the in- 
completeness fraction are at the level of 1 per cent and the 
error is dominated by the systematic uncertainty associated 
with the dc offset (Section 4.3.1). The systematic uncer- 
tainty becomes less significant for the EW> 2.0 A sample. 

Nestor et al. (2005) find the the EW distribution of a 
sample of 1 300 absorbers (with redshifts 0.366 < z < 2.269) 
from SDSS DR4 is well described by an exponential dis- 
tribution. They find that although moderately strong lines 
(0.4 A < Wq 2796 < 2.0 A) show no evidence for redshift 
evolution, the absorbers with larger EWs show an increase 
in number with redshift that is more pronounced for ab- 
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Figure 15. The flux ratio against wavelength for a composite mean spectrum created from 18 objects satisfying Equation 10. The 
absorbers have been chosen (using the D m ax statistic) as those most likely to exhibit LMC-type dust. The underlying composite spectrum 
is best fit by the LMC composite extinction curve, with mean E(B — V)smc = 0.16, E(B — V)lmc — 0-21 and E(B — V)mw — 0.18. 
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Figure 16. The flux ratio against wavelength for a composite median spectrum created from 1 792 absorbers with ^ E(B — V)smc ^ 
0.05 and z a b s > 0.9. The shape of the flux ratio spectrum for the low E(B — V) composite is consistent with SMC-type extinction 
with median E(B — V)smc = 0.018, E(B — V)lmc — 0.022 and E(B — V)gal = 0.016. The minimum flux ratio on the y-axis has 
deliberately been set to 0.84 and the x-axis extended to 4100 A to better show the shape of the extinction curve at small deviations from 
unit flux-ratio. A composite based on the mean flux-ratio for the sample is noisier but extremely similar in shape. 



sorbers of higher EW. Nestor et al. interpret this redshift 
dependence as an evolution in the kinematic properties of 
absorbers over a range of intermediate redshifts. Lundgren 
et al. (2009) find a similar result which suggests stronger 
evolution in the redshift number density of strong Mgll ab- 
sorption systems relative to lower EW samples, albeit over 
a smaller redshift range 0.4 ^ z ^ 0.8. 

Using our parametrized dust correction (equation 7), 
we predict that the observed relative fraction of high EW 
systems will decrease by 36 per cent over the redshift inter- 
val 0.4 < z < 2.0. Our results thus predict that the true 
increase in the proportion of high EW systems with increas- 
ing redshift is significantly greater than found by Nestor 
et al. (2005) from the observed population of absorbers. 



5.5 The GRB/Quasar absorber number density 
sight line discrepancy 

We have discussed the effects of dusty high EW absorbers on 
the completeness of optical magnitude-limited quasar sam- 
ples over a range of absorber redshifts. The highest red- 
shift absorbers have been observed using GRB optical after- 
glows as background sources (Olivares et al. 2009). Although 
GRBs would individually experience the same obscuration 
effects, GRBs are found to have a huge range in intrinsic 
brightness and spectra are obtained at variable time inter- 
vals following the GRB's peak brightness. The consequent 
very large dispersion in apparent GRB brightness, rather 
than the presence of modest amounts of dust in any inter- 
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Table 1. A sample of 25 objects with E(B - V) S mc > 0.3 



SDSS name 


•Zqso 


^abs 


WA2796 

(A) 


E(B - 


^)SMC 


E(B- 


- V) F itz 




Yv e ga 


J\ega 


Hvega 


KVega 


J 081524. 62+555153. 3 


1 


781 


1.32d(l.o08) 


1.9(1.3) 





37 





31 


18.63 










J083257. 64+333214. 6 


1 


017 


0.716 


2.2 





33 





29 


19.06 










Tnoro a a >t a i o /for" a r\ a n 

J 085244. 74+343540.4" 


1 


656 


1.310(0.884) 


3.1(2.5) 





31 





27 


18.21 










J091753. 89+000300. 8 


2 


139 


0.729 


4.3 





34 





31 


18.78 


17.92 


17.46 


16.69 


15.84 


J 092053. 12+385002.7 


1 


276 


0.472 


1.8 





30 





28 


17.96 










J 092333. 29+595747. 8 





839 


0.662 


2.6 





31 





28 


19.02 










J 093444. 30+172453. 5 


1 


899 


0.553 


2.8 





33 





32 


19.00 










J 093738. 04+562838. 9" 


1 


804 


0.978 


4.9 





30 





23 


18.50 










J100713. 68+285348. 4 C 


1 


050 


0.884 


3.5 





60 





50 


18.25 










J 113152. 15+435318.3 


2 


139 


1.098(1.377) 


2.4(1.1) 





34 





28 


18.82 










J 113811.59+382119.5 





899 


0.759 


1.9 





44 





38 


18.41 










T 1 oaoai o o i c\r> a a a v n, 

J 120301.00+063441.5" 


2 


180 


0.862 


5.6 





42 





37 


18.44 


17.28 


16.98 


16.41 


15.51 


t I ooo i o r* a i /f oor\or\ o 

J 120913.61+433920.9 


1 


397 


0.412 


2.4 





30 





28 


17.35 










J 121547. 11+293409.9 


1 


532 


0.788 


1.0 





30 





25 


17.88 










t -i c\ a r\ a r* o i -i oo a r\r\ r\ 

J 124946.59+120400.0 


1 


688 


1.469 


2.3 





33 





31 


18.20 


16.90 


16.28 


15.49 


15.07 


J131103.19+551354.3 





927 


0.600 


1.2 





39 





37 


18.87 










J133753.11+052746.2 


1 


743 


0.575 


1.4 





31 





30 


17.99 






16.05 


15.72 


J140807.05+374457.3 


1 


497 


0.776 


3.0 





30 





25 


18.61 










J144611. 68+484613.7 


1 


202 


0.669 


2.2 





38 





34 


18.46 










J144642.92+012552.4 


1 


422 


0.525(0.510) 


2.4(1.4) 





34 





32 


18.08 


17.48 


16.92 


16.13 


15.71 


J145344.23+102557.5 


1 


770 


0.757 


2.6 





38 





34 


19.02 


18.09 


17.66 


17.10 


16.40 


J145907.19+002401.2 6 


3 


037 


1.394 


1.7 





31 





23 


17.83 










J155435. 08+484411.1 


2 


041 


0.912(1.385) 


2.4(1.1) 





31 





24 


18.98 










J170220.06+591538.6 


1 


798 


0.724 


1.8 





32 





28 


18.73 










J171123.03+311613.7 


2 


041 


1.256(1.697) 


2.3(1.9) 





31 





27 


18.65 











Magnitudes are corrected for Galactic extinction using the prescription of Schlegel et al. (1998). Photometric errors for the SDSS i-band photometry 
are normally < 0.03 mag and for the UKIDSS photometry ~0.05 mag. The absorber redshifts and EWs in brackets correspond to those spectra which 
have two Mgll systems with EW>1 A. Errors in the rest-frame EWs are typically ±0.1 A. E(B — V)pitz is the result of the fit to the extinction 
curve of Fitzpatrick (1999) with Ry =2.1 (see Appendix A). Systems that have been the subject of previous investigations are designated ( a ) 
Srianand et al. (2008), ( b ) Wild & Hewett (2005a), ( c ) Zhou et al. (2010), ( d ) Wild et al. (2006) and ( e ) Jiang et al. (2010). 
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Figure 17. Left-hand panel: The EW distribution of the absorber sample (solid-line), and corresponding dust-corrected EW distribution 
(dashed-line), which has been calculated using equation 6. Right-hand panel: The redshift distribution of the absorber base sample (solid- 
line), and the corresponding dust-corrected distribution (dashed-line). 
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vening absorbers, dictates whether an object is observed and 
absorption systems identified. Therefore, one would expect 
a larger number of moderate E(B — V) absorption systems 
to be present in GRB spectra than quasar spectra. 

Prochter et al. (2006) identified 14 strong intervening 
Mgll systems (at a mean redshift of (z) = 1.1) along 14 
GRB sight lines, an incidence roughly four times higher than 
along sight lines to quasars. The result is not expected if 
both GRBs and quasars sample random lines of sight. Since 
the intervening absorption systems are thought to be inde- 
pendent of the background source, the observed discrepancy 
has led to a call to review the fundamental assumptions that 
underpin extragalactic absorption line research. The discrep- 
ancy was confirmed by Sudilovsky et al. (2007), but the am- 
plitude has since been reduced to a factor of 2.1 ± 0.6 using 
a larger sample of 22 absorbers (Vergani et al. 2009). 

A series of explanations have been suggested to explain 
the observed GRB/quasar discrepancy. Porciani et al. (2007) 
claim that strong Mg II absorbing gas may be intrinsic to the 
GRB circumburst environment or originate from supernova 
remnants lying in the same star-forming region. It is also 
claimed (Prochter et al. 2006; Porciani et al. 2007; Tejos 
et al. 2009) that source amplification due to strong gravi- 
tational lensing may bias the GRB spectral samples toward 
targets that contain more intervening absorbers. Another 
mechanism has been proposed (Frank et al. 2007), which 
claims that the discrepancy is due to the different beam sizes 
of GRBs and quasars, but this has subsequently been ruled 
out by observational analysis (Pontzen et al. 2007). Finally, 
a number of authors (Prochter et al. 2006; Porciani et al. 
2007) have suggested that dust associated with strong Mg II 
systems results in a reduction in the observed number of ab- 
sorbers seen along sight lines to quasars in flux-limited cat- 
alogues. In general, most authors agree that the differences 
between Mgll toward quasar and GRB sight lines cannot 
be due to a single effect, but dusty absorbers have not been 
thought to be a significant factor. 

In light of our findings on the existence of a signifi- 
cant population of dusty Mgll absorbers, we can provide a 
more definitive determination of the effect of dust on the 
GRB / quasar sight line discrepancy. Using a sample of GRB 
sight lines with intervening absorbers we can assess how this 
sample would suffer from incompleteness if the sight lines 
were instead illuminated by quasars. We assume that the 
GRB sample is unbiased and probes all intervening Mgll 
absorbers up to a moderate absorber E{B — V) of 0.2 mag. 
Given an EW and redshift distribution of Mgll absorbers 
in GRBs similar to our quasar sample (Fig. 17), we would 
expect the GRB sight lines to show the missing 24 per cent 
of absorbers (Section 5.4). However, in general GRB spectra 
contain absorbers with different redshift and EW distribu- 
tion to the absorbers in quasar spectra. The observed dis- 
crepancy in number density due to dust will be sensitive to 
such differences. 

Two GRB absorber samples in the literature are the 
Vergani et al. (2009) sample and the Fynbo et al. (2009) 
sample. The Vergani et al. sample contains 22 strong systems 
(EW> lA) identified using high-resolution UVES sample 
and a series of high- and low-resolution GRB afterglow spec- 
tra from the literature. The Fynbo et al. sample contains 15 
strong systems present in the follow-up spectra of 77 op- 
tical afterglows of Swift detected GRBs. The redshift dis- 



tributions of the two samples are shown in Fig. 18-(left). 
The redshift distributions for the two GRB samples both 
exhibit a larger fraction of absorbers at high redshift com- 
pared with the quasar absorber redshift distribution. Fig. 
lS-(right) shows the degree to which both samples suffer 
from EW incompleteness by comparing the EW distribu- 
tions with an exponential (Nestor et al. 2005) of the form 

N (W X2796 ) oc e -< 2796 /w (11) 

where W* is the exponential scale factor. A fit to the cor- 
rected quasar absorber EW distribution (Fig. 17) gives a 
value of iy*=0.746±0.008. Scaling the corrected quasar EW 
distribution using the redshift paths appropriate to the two 
GRB-derived samples results in the predicted number ver- 
sus EW distributions for the Vergani et al. and Fynbo et al. 
GRB samples shown in Fig. 18-(right). Both samples are 
expected to be incomplete at low EW but the scaled dis- 
tribution functions accurately describe the Vergani et al. 
sample with EW>2.0 A, and the Fynbo et al. sample is con- 
sistent over the interval ~ 2. 0-5.0 A. Our flux- limited quasar- 
derived Mgll absorber distribution has no information re- 
garding high-EW systems with Wq A2796 >5.0 A (where the as- 
sociated large E(B — V) essentially removes nearly all such 
systems from the SDSS quasar sample). The extension of 
the distribution of Mgll absorbers beyond Wq 2796 =5.0 out 
to ~10A is clear and strongly suggests that the absorber 
dust-induced differences between quasar- derived and GRB- 
derived samples will be even larger than calculated using 
our Af(Wo A2796 ) parametrization. 

We can calculate the fraction of absorbers that would be 
missed in corresponding flux-limited quasar sample by esti- 
mating an E(B — V) according to equation 7 and assigning a 
weighting according to the dust obscuration bias/signal-to- 
noise ratio effects (equation 6). The predicted discrepancy 
is obtained by calculating the fraction of absorbers that are 
missed from a quasar- derived sample with an Wq 2796 dis- 
tribution corresponding to Equation 11 (with W* — 0.746) 
and redshift distributions of the Vergani et al. and Fynbo 
et al. (2009) absorbers 16 . 

The strong dependence of the dust content on the EW 
of the Mgll absorbers means that the degree to which the 
incidence of absorbers is affected is a very strong function 
of the low EW limit of the samples (where the majority of 
absorbers reside). Fig. 19 shows the predicted discrepancies 
for a sample with the same absorber redshift distributions as 
Vergani et al. (2009) and Fynbo et al. (2009) as a function of 
the lower EW limit for the absorber sample. The increasing 
impact of dust, as samples with increasingly large lower EW 
limits are considered, illustrates the importance of incorpo- 
rating quantitative estimates when inter-comparing results 
from different observed samples. The discrepancy reaches a 
full factor of two for W^ 2796 >2.4 A. 

Our results are consistent with the investigation of Tejos 
et al. who find no significant increase in number density 
of low-EW absorbers, EW<1 A, towards GRBs. For GRB 
samples where high-EW absorbers (EW>2A) represent a 



16 Note that the procedure adopted is inherently conservative 
because the parametrized dependence of the E(B — V) values on 
Wq 2796 reduces the impact of individual dusty absorbers. 
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significant fraction of the absorber sample, our predicted en- 
hancement in Mg II number density is a factor two compared 
to absorber samples seen towards SDSS-selected quasars. 
The evidence from the Fynbo et al. sample for the exis- 
tence of a significant number of absorbers with W A2796 s as 
high as 10 A suggests that the true discrepancy will be sig- 
nificantly greater than shown in Fig. 19 for samples with 
minimum detectable Wo 2796 above 1 A. The amplitude of 
the enhancement is higher than has been reported in the 
literature hitherto and dusty absorbers could be responsible 
for a significant proportion of the absorber number density 
enhancement seen towards GRBs. 

A careful study of the incidence of Mg II absorbers to- 
wards a third class of objects, the BL Lacs (or Blazars), 
whose selection is closely akin to GRBs in not being flux- 
limited in the optical, has been conducted by Bergeron et al. 
(2011). They find an excess of Mgll absorbers along Blazar 
sight lines compared with quasar sight lines of ~2.0 with no 
apparent difference between low-EW absorbers (EW<lA), 
and high-EW absorbers (EW>1 A). However, the absorber 
sample is numerically small and while strong differences be- 
tween low- and high-EW absorbers are not present, less sig- 
nificant differential effects as a function of absorber EW are 
not ruled out with any confidence. 

Our estimate of the true number density of strong 
(EW>1 A) Mgll absorbers along sight lines to SDSS quasars 
is some 15 per cent higher than the value used by Berg- 
eron et al. (2011), reducing the observed overdensity some- 
what. The median redshift of the Blazar absorber sample is 
low, even lower than for the SDSS quasars sample presented 
here (z=0.84 compared to £=1.1), and the effect of dusty 
absorbers on the statistics is correspondingly much lower 
than for the GRB absorber samples discussed above. Our 
model predicts an increased overdensity for the Bergeron 
et al. (2011) sample of only 15 per cent. We also note that 
the contribution to the overdensity signal may be redshift 
dependent (fig. 2 of Bergeron et al.) with the modest path- 
density at z>1.0 producing much of the effect. The Bergeron 
et al. (2011) study does not offer direct support in favour 
of the predictions from our dusty absorber model but nor 
does it rule out such a contribution. However, assuming that 
the factor of two overdensity of absorbers at low redshifts 
for Blazar samples is confirmed, then more contributors, in 
addition to the effects of dusty absorbers, are required to 
explain the observed absorber overdensity towards GRBs. 



6 DISCUSSION AND CONCLUSIONS 

With the aid of the enormous statistical power of the SDSS, 
we have been able to provide an improved characterisation 
of the dust content of Mg II absorption systems and its ef- 
fects on the properties of absorber samples obtained using 
magnitude-limited quasar surveys. The principal conclusions 
of our investigations are: 

(i) We find that dust content (E(B — V)) in a large sam- 
ple of ~8300 strong Mgll absorbers increases much more 
strongly as a function of absorber rest-frame equivalent 
width (EW) (see equation 7) than previously reported. 

(ii) The effects of dust on the observed Mgll absorber 
population results in significant changes to the statistics of 
the intrinsic absorber population. We find 24 ± 4 per cent of 




^ I i i i i i i i i i i i 

0.5 1 1.5 2 2.5 3 
lower EW limit 



Figure 19. The predicted discrepancy between GRB and flux- 
limited quasar sample sight lines due to dust as a function of lower 
Wq A2796 limit, for a GRB sample with EWs distributed accord- 
ing to Equation 11 (with W* = 0.746), and redshift distribution 
consistent with that of the Vergani et al. (2009) and Fynbo et al. 
(2009) samples. 

absorbers are missing from the sample for EW > 1.0 A and 
34 ± 2 per cent for EW > 2.0 A. For absorbers with EW 
> 3.0 A) the fraction rises to more than a factor of two. 

(iii) We find no significant detection of an evolution in 
absorber dust content over the redshift interval 0.4^2^2.2 
but the constraint is relatively weak, with a la uncertainty 
of 0.3 in the power-law index a (E(B - V)oc(l + z) a ) 

(iv) The form of the extinction curve for absorbers is 
found to be of SMC-type at low values of E(B — V). For 
absorbers with modest E(B — V)<0.2 at redshifts z a bs > 1, 
where the 2175 A is visible in the SDSS spectra, approx- 
imately a third of the population shows evidence for the 
presence of dust with an LMC-like extinction curve. 

(v) For absorbers with the highest E(B — V)>0.3 there is 
evidence for an extinction curve that differs from the SMC, 
LMC and MW extinction curves. The high E(B— V) systems 
possess a 2175 A feature similar in strength to that seen in 
the MW but the overall shape of the extinction curve in the 
near-ultraviolet through to the near- infrared, is greyer than 
any of the SMC, LMC or MW extinction curves. 

(vi) The impact of dust on the observed discrepancy in 
the Mgll absorber sight line densities between GRBs and 
quasars depends on both the absorber redshift distribution 
and the minimum equivalent width limit of the samples ob- 
tained towards GRBs. Samples with high lower EW-limits 
for the detection of Mg II absorbers and significant redshift 
path for absorbers at high redshifts are most affected. Pub- 
lished GRB samples are predicted to possess absorber den- 
sities a full factor of two greater than for quasars from flux- 
limited samples due to the effects of dust extinction. 

Higher EW Mg II systems are found to exhibit a higher 
star formation rate (Menard et al. 2009) and the dependence 
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Figure 18. Left-hand panel: The redshift distribution of the Vergani et al. GRB absorber sample (solid-line), and the Fynbo et al. 
sample (dashed-line). The absorbers exhibit similar redshift coverage, although the Fynbo et al. sample shows a slightly larger fraction 
of absorbers at higher redshift than the Vergani et al. sample. Right-hand panel: The EW distribution of the Vergani et al. (2009) GRB 
absorber sample (solid-line), and the Fynbo et al. (2009) sample (dashed-line). Overlayed are the results of an exponential fit (Equation 
11 with W* = 0.746) to the two EW distributions. The vertical scaling is calculated using the redshift path of the two samples. Both 
samples are expected to be highly incomplete at low EW. The extension of the absorber distribution to Wq 2796 ~10 A in the Fynbo et al. 
sample is clear. 



of E(B — V) on EW can be explained in terms of the star for- 
mation rate of the absorbing galaxy. Supernovae associated 
with star formation impart mechanical energy to the cold 
Mg II absorbing gas, which leads to velocity spreading of the 
absorption feature and a consequent increase in the observed 
absorber EW. An increased amount of star formation also 
leads to increased dust production. The presence of a larger 
amount of dust in absorbers associated with the highest star 
formation rate has potential consequences for determining 
the redshift-evolution of absorbers with the highest EWs al- 
though the impact on the overall star formation rate density 
of the Universe (Hopkins 2004) is unlikely to be more than 
a perturbation. 

Explaining the observed differences between the Mgll 
absorber density towards quasars in flux-limited samples, 
GRBs and BL Lacs requires considerably larger absorber 
samples to map out the dependence of the differences on ab- 
sorber EW and redshift. While the explanation will almost 
certainly involve more than one effect, the results presented 
here show that obscuration due to dust associated with the 
absorbers is almost certainly a significant factor. 

To date, the limited statistics available has meant that 
constraints on the dust content and form of the extinction 
curves associated with metal absorbers have been weak. The 
availability of metal absorber samples with ^10 000 sys- 
tems and the detection of a small number of systems with 
E(B — y)>0.3 is changing the situation, as shown in this 
paper. At low values, E(B — y)<0.05, the extinction curve 
appears to be essentially featureless in the near- ultraviolet 
with an SMC-like form. A new result from our study is 
the detection of a weak 2175 A feature in a significant frac- 
tion of Mgll absorbers with E(B — V)~0. 1-0.2, consistent 
with an LMC-like extinction curve. Individual systems with 
higher E(B — V) values of >0.3 and an apparent strong 



2175 A feature have been identified in recent years. Our re- 
sults for the Mgll absorbers with the highest E(B — V) 
values, particularly the small subset with near-infrared pho- 
tometry available, indicate that there may be a systematic 
dependence in the form of the extinction curve as a func- 
tion of E(B — V) and, more fundamentally, the column den- 
sity of the absorbers. Obtaining near-infrared photometry, 
to combine with SDSS spectra, for a larger sample of such 
objects is necessary to establish the general nature of any 
such trend. Confirmation of a systematic change in the form 
of the extinction curve as advocated here, from featureless 
and SMC-like at low E(B — V), through to greyer overall 
than any of the SMC, LMC or MW curves but with a strong 
2175 A feature (consistent with the low Rv form presented 
by Fitzpatrick (1999)) at high E(B - V), would provide 
powerful constraints on models of star formation and chemi- 
cal enrichment histories for physical systems associated with 
metal absorbers over an extended range in column density. 
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Table Al. The Fitzpatrick (1999) extinction curve with R v =2.1, 
over the wavelength range 1000 - 25 000 A. 
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Figure Al. The Fitzpatrick (1999) extinction curve with Ry = 
2.1. 
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